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Abstract 

In non-coherent wideband fading channels where energy rather than spectrum is the limiting 
resource, peaky and non-peaky signaling schemes have long been considered species apart, as the 
first approaches asymptotically the capacity of a wideband AWGN channel with the same average SNR, 
whereas the second reaches a peak rate at some finite critical bandwidth and then falls to zero as 
bandwidth grows to infinity. In this paper it is shown that this distinction is in fact an artifact of the 
limited attention paid in the past to the product between the bandwidth and the fraction of time it 
is in use. This fundamental quantity, called bandwidth occupancy, measures average bandwidth usage 
over time. For all signaling schemes with the same bandwidth occupancy, achievable rates approach to 
the wideband AWGN capacity within the same gap as the bandwidth occupancy approaches its critical 
value, and decrease to zero as the occupancy goes to infinity. This unified analysis produces quantitative 
closed-form expressions for the ideal bandwidth occupancy, recovers the existing capacity results for 
(non-)peaky signaling schemes, and unveils a trade-off between the accuracy of approximating capacity 
with a generalized Taylor polynomial and the accuracy with which the optimal bandwidth occupancy 
can be bounded. 
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Index Terms 

Wideband regime, peaky signals, non-peaky signals, bandwidth occupancy 

I. Introduction and Related Work 

Recently there has been great interest in wireless channels with a large bandwidth, owing 
in part to the prospective investments onto the millimeter wave bands, where vast quantities of 
new spectrum is readily available [|2l-[[5l|. In a frequency selective fading channel where there is 
no channel state information at the receiver (CSIR) or the transmitter, the wideband capacity is 
affected by the growing uncertainty in the channel impulse response. As bandwidth grows while 
energy is constrained, it becomes infeasible to estimate the channel coefficients to a precision 
sufficient for coherent detection. Moreover, if the transmitted signal power is spread across all the 
available bandwidth and time slots, the desired signal would be buried by the channel uncertainty 
when bandwidth is too large. Medard and Gallager proved this [O through an upper bound to 
the rate that is proportional to the ratio between the fourth moment of the signal (E [|x|^]) and 
its bandwidth {B), i.e., R <oc E [|x|^]/i?. That is, to achieve rates above zero when B ^ oo, 
one has to make E [Ixl"^] grow at least as fast as B by concentrating the power of the signal in 
a vanishing fraction of its transmitted symbols (i.e. infrequent bursts of very large power). 

In this paper we investigate the capacity bounds of non-coherent wideband fading channels in 
multi-input multi-output (MIMO) setup where both the signaling bandwidth and signal peakiness 
are design parameters. The channel is assumed to be rich scattering, frequency selective, block 
fading with a coherence time Tc and a delay spread D, such that the channel frequency response 
becomes uncorrelated for frequencies apart from more than one coherence bandwidth Bc^l/D. 
The channel coherence length, BcTc, is assumed to be large for capacity analysis purposes, as 
in almost all practical channels, BcTc 3> 1. In our expressions we temporarily treat BcTc as a 
fixed parameter to derive closed-form expressions, where approximation errors originated from 
are highlighted in small-o expressions parametrized by higher order terms of BcTc. 
We further assume that BcTc>Nt, which is easily satisfied in typical systems where the number 
of transmit antennas is not massive. We generalize the analysis method in [|71, developed for 
non-peaky signaling in single-input single-output (SISO) systems, to MIMO systems and extend 
it to arbitrary level of signal peakiness by enforcing a transmission duty cycle 5 g( 0, 1]. The duty 
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cycle prescribes a bursty transmission seheme where the transmitter is aetive only for a fraetion 
S of time with boosted signal power P/S harnessed from the (1—5) silent-eyele. Denoting by 
C(B) the eapaeity of the uneonstrained non-eoherent ehannel and by C(B, S) the maximal rate 
aehieved by using bandwidth B and duty eyele S, for all S > 0 and S G (0,1], we have 

C(B,S) < C(B) < = NrP/No [nats/s], 


where C'°° is the limit eapaeity of the eoherent ehannel at infinite bandwidth, P is the reeeived 
signal power, No is the noise power speetral density, and Nj. is the number of reeeive antennas. 
Note that the first inequality is striet beeause we do not exploit the position of the aetive symbols 
to eonvey information. We show in See. [Ill] that C(B, S) is upper and lower bounded by 

R^^(SB) < C(B,S) < R^^(SB). 


Note that both the upper and lower bounds, up to a small approximation error o(l/SB), depend 
on B and S only through the produet SB, whieh measures average bandwidth usage over time 
and is named the “bandwidth oeeupaney”. Our results show that for a series of signaling sehemes 
with finite signaling bandwidth B larger or equal to a critical bandwidth occupancy {SB)crit, 
whieh falls in a range preseribed by elosed-form expressions, it is possible to aehieve rates elose 
to within the same rate penalty 


C{B,S)>K—-Ac, 

i\o 

as long as the duty eyele is 5= 


Ac = nPJ '-^ ajf (k - 2 + iV. + iV,) log ,r, ( 1 ) 


. (SB), 


^ Here iVt is the number of transmit antennas and fi:>0 
is the kurtosis (whose definition is deferred to See. HH) of the ehannel. Thus, it is possible to 
approaeh C°° up to the same gap with any (5g(0, 1]. Note also that BcTc:>l leads to Ac^O and 
Furthermore, we show in See. UV] that the analysis of C{B) with peaky signaling 
in literature |[8l, flUl experienees exaetly this same gap to C°°, although we obtained Q using 
non-peaky signals 0 and a power-boosting duty eyele (5e(0,1]. Fig. [U illustrates the relation 
between our bounds C{B,S), eapaeity C{B), and the eoherent wideband ehannel limit C'°°. 

The main eontribution of this paper is the unified approximation of with peaky and non- 
peaky signaling, showing that these two extremes ean be eonneeted by all level of peakiness 
parametrized by the duty eyele 5g(0, 1]. All signaling sehemes {B, 5) with the same bandwidth 
oeeupaney SB={SB)crit approaeh within the same eapaeity gap up to a small approximation 
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Figure 1. All transmission strategies with the same bandwidth occupancy SB={SB)ciit achieve the same polynomial 
approximation of C°° at different bandwidths. C{B) is separated from the maximum C{B, S) by a difference of o(l/((5i3)crit)- 


error of o(l/((5i?)crit)- We have also derived elosed-from expressions for eapaeity bounds and 
eritieal bandwidth oeeupaney for all values (5g(0, 1], whieh provide valuable engineering insights 
and tools to quantify the resources needed to approach C°°. As a byproduct, we obtained a group 
of closed-form bounds to the range of {6B)crit that are implicit in the existing literature [[8]|, ll9ll . 
These parametric bounds can be tuned based on an accuracy-resolution tradeoff to complement 
the range identified in our non-peaky signaling analysis. 

A. Related Work 

The results in [|6l have been extended to signals with output fourth-order constraint [fTOl or 
small input peakiness constraint [QTI . Telatar and Tse [fT^ related channel uncertainty to the 
number of resolvable independent paths, and showed that in a rich scattering environment where 
this number increases with B, the rate can grow as long as the signal power in each path is not 
too low, but it starts decreasing when the number of paths is above a critical value. 

The capacity of a wideband fading channel achieves first order optimality if, as B goes to 
infinity, it has the same limit as a wideband additive white Gaussian noise (AWGN) channel. 
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This has been studied in [[8]|, JHl, ITT^ - IIT?! and the linear in power capacity limit for MIMO is 


lim C{B) 

B^oo 


noncoherent 


= lim C{B) 

B^oo 


AWGN 


N,P 

1^' 


To quantify the “exchange rate” of bandwidth to capacity in the asymptotic regime where 
B —)■ oo, the concept of wideband slope was introduced in [[T3l . A larger wideband slope means 
that higher rate gain is obtained given the same amount of extra bandwidth. The wideband slope 
is studied in [fT3l Theorem 9] based on the second order term of a Taylor series expansion of the 
spectral efficiency {C/B, in nats/s/Hz) with respect to the signal-to-noise ratio (SNR) at each 
receive antenna, SNR=P/(i?A"o). The wideband slope is inversely proportional to the second 
order derivative of the spectral efficiency at SNR=0, which is finite for AWGN and coherent 
fading channels (i.e., with perfect CSIR) but —oo for non-coherent scenarios (i.e., with no CSIR). 
Thus the coherent fading channel is second order optimal but the non-coherent channel is not. 
This abrupt distinction contrasts with the intuition that, as the channel coherence time T^. and/or 
frequency Be grow, channel estimation becomes increasingly rewarding and the capacity of 
the non-coherent channel converges to the capacity of the coherent channel. This contradiction 
was resolved in flSl, BH by showing that in non-coherent Rayleigh fading channels the spectral 
efficiency C{B)/B is better represented by a generalized Taylor polynomial of order 1-l-a < 2, 


C{B) 

B 


--K 


P 


N,{N,+Nt) f P 


2N, 


NnB 


l-\-a 


+ o( 


51+c 


), [nats/s/Hz], 


( 2 ) 


where the exponent Q!G(0, 1) grows with increasing BcTc. The first term equals C°°/B, repre¬ 
senting a first order optimal upper bound of the spectral efficiency when rate is power-limited. 
The third term captures the approximation error, that vanishes faster than as B^oo. 

The second term represents the penalty from lack of channel knowledge. It contains SNR^"*"", a 
sub-quadratic term (l-|-a<2) that characterizes the convergence speed of the spectral efficiency 
for non-coherent fading channels. Representing Q by the second order Taylor polynomial leads 
to an infinite coefficient to the term SNR^ (wideband slope) and lack of second order optimality 
as in lfT3l . In this paper the word “polynomial” refers to these generalized Taylor polynomials 
with real-valued exponents. 

Although peaky signals are imperative to achieve first order optimality ifT^ Th. 7], they are 
challenging to synthesize owing to hardware non-linearity and the infinite amount of bandwidth 
they require in non-coherent channels. If a small gap from at a large but finite bandwidth is 
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admissible, which is the case in all practical applications of asymptotic results, recent works have 
shown that non-peaky signals may suffice. For example, Zhang and Laneman ifTSll investigated 
the achievable rate of phase-shift keying (PSK) for frequency-flat time-varying non-coherent 
Rayleigh fading channels. Under average power constraints, this signaling scheme approaches 
the wideband capacity limit for low but not too low SNRs. For signals subject to both peak 
and average power constraints, it was observed in ifT^ that the gap between capacity upper and 
lower bounds can be very small for discrete-time frequency-flat Rayleigh fading channels. The 
capacity of non-coherent time-frequency selective wide-sense stationary uncorrelated scattering 
(WSSUS) channels with both peak and average power constraints has been studied in [fTTIl . 
where bell-shaped capacity upper and lower bounds were established and the capacity optimal 
bandwidth, the critical bandwidth, was coarsely identified as a function of the peak power and 
the scattering function. For flat scattering functions, the capacity bounds depend on the system 
bandwidth and the input-signal peak constraint only through their ratio. The results in ifTTl have 
been extended to MIMO in [|T^ . where the impact of transmitter/receiver antenna correlation on 
capacity was also investigated. Lozano and Porrat [[71 considered non-peaky signaling in SISO 
systems under a general fading distribution. Their results show that, when bandwidth is not 
too large, there is a transitory first stage where rate R{B) grows with B before approaching a 
maximum i?(i?crit) at the critical bandwidth B=Bcn^, beyond which the rate R{B) decreases to 
zero as B grows unbounded. By resorting to computation of mutual information rather than the 
capacity analysis as in [fTTl . [fT^ . they provided closed-form expressions to the maximum rate 
and the corresponding capacity gap. 



(3) 


where A vanishes with increasing coherence length B^Tc- For Rayleigh fading, closed-form 
expressions for the range of Bent were also derived. 

Even though [fT3l Th. 7] found that peaky signaling is imperative to achieve first order 
optimality, the definition of first order optimality enforces an implicit requirement to make 
bandwidth grow as high as possible (B —)■ cxd). Thus, only those inputs that approach C°° when 
B is infinite are covered by lfT3l Th. 7]. What our results show is that C°° can be approached 
as well using a finite bandwidth B and non-peaky signaling. 

Unlike in [j8l, dH where the non-coherent wideband fading channel capacity C{B) is obtained 
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by using the position of signal pulses in the frequency domain (i.e., FSK) to convey information, 
in our analysis the position of actively transmitted symbols in the time domain, which collectively 
defines the active-cycle, is revealed in advance to the receiver and therefore bears no information. 

Our capacity bounds are based on computation of mutual information with constrained input 
signal peakiness - in the sense of kurtosis - that is controlled by enforcing a duty cycle 6 E (0,1]. 
This is in contrast to lUTl . lITSl where capacity analysis is used with peak constraint on the 
amplitude of transmitting signals. Our choice of mutual information analysis can be justified 
from two aspects: even thought we do not design inputs to achieve the capacity bounds we 
can guarantee such inputs exist as long as the channel and noise are stationary weakly mixing 
processes, see lfT9l Prop. 2.1]; the rate upper and lower bounds and the range of the critical 
bandwidth occupancy can be described in closed-form expressions, which are otherwise difficult 
to obtain using capacity analysis, see ifTTIl . ifT^ . 

Our choice of using duty cycle rather than peak constraint on signal amplitude [|T^ . ifTTl to 
control the signal peakiness can be justified as follows: given the same average power constraint, 
a peak constraint on signal amplitude will limit the peak-to-average power ratio (PAPR), which 
is sufficient but not necessary to generate a constraint on signal peakiness. Signals with finite 
peakiness may have infinite PAPR (e.g., Gaussian signal has infinite PAPR but only a small 
kurtosis k=2). It must be noted that in non-coherent wideband fading channels, capacity is 
related to the peakiness in the kurtosis sense [1^, IHAII . 

The rest of this paper is organized as follows. We introduce the system model in Sec. HI] and 
present our unified analysis of wideband non-coherent channel in Sec. Hill We describe our non¬ 
coherent polynomial approximation to coherent capacity, and discuss its relation with literature 
in Sec. HYl Finally our conclusions are in Sec. jV] 

II. System Model 

We consider a rich scattering, frequency selective, block fading. At x Nr MIMO wideband 
channel with an impulse response between antennas (m, v). For compactness we assume 

that all channels experience a coherence time T^. and a delay spread D and the channel fre¬ 
quency response becomes uncorrelated for frequencies apart more than one coherence bandwidth 
Bc^l/D. We focus only on the frequency signaling scheme since it is known [|71 that differences 
between frequency and time signaling only affect the scaling with bandwidth in its vanishing 


higher order terms. In the following we present the eharacteristics of the discrete-time system 
modeli]. Justification of our choice of the wideband fading model is presented in Appendix 
Our model starts from a continuous-time wideband fading channel, followed by the discretiza¬ 
tion/sampling process on the input-output signals. This provides a signaling scheme where every 
Tc seconds, the transmitted signal x^'^'^[n] with bandwidth B carries K=BTc complex samples 
on antenna mg[ 0, At—!]• Taking a A-point DFT of the complex samples for each antenna and 
then stacking all the At vectors up, the transmitted codeword is uniquely defined by the At A x 1 
vector X that satisfies the average power constraint 

1 e [1x1=] < Pn. 


For i=kNt+u, the i-\h coefficient of x, denoted as x^'^\ corresponds to the transmitted signal 
on antenna u with DFT index fcG{0,1, ..., A— 1}. For each pair of antennas {u,v), the dis¬ 
crete samples of the channel have M=BD i.i.d. coefficients n=0,1, ... ,M— 1, with 

M/K=D/Tc=^^. After applying A-point DFT to each discrete channel sequence 
we define a block-diagonal matrix 

/ wrm n ... n \ 


H 


H[0] 

0 


0 

0 

H[l] 






0 

0 


0 

H[A - 1] 


(4) 


where H[/c] contains in its (u, M)-th element the fc-th DFT coefficient of [n], whose distribu¬ 
tion is determined by the impulse response Each channel only has M i.i.d. coefficients 

and any two blocks H[fc] and H[/c'] are correlated only if \k—k'\<BcTc. We also define the 
average gain of the n-th channel coefficient gY‘’'^'^=E [n] p] satisfying 9^’^^ = 1- 

When D<tiTc, a cyclic prefix with negligible influence in rate can be inserted to remove the 
inter-symbol interference and the signal received on each fading realization, Tc, depends only 
on the state of the channel and signal transmitted during the same realization. After applying 
A-point DFT to the received signal, we can represent the system as 


y = Hx -f z, 


(5) 


*The equivalence between the discrete-time and continuous-time channel models for SISO is established in (201 using sampling 
and DFT, and in (H] using pulse shaping filter banks with Weyl-Heisenberg projection. Our result uses MIMO in a rich scattering 
environment and we provide explicit mapping of the channel coefficients between two different discrete-time models. 











where y is a N,.K x 1 vector whose i-th element with i = kNj.+v, corresponds to the signal 
received on antenna v with DFT coefficient index k. The noise vector z follows a Gaussian 
distribution with PSD Nq (CAA(0, In^kNqTc)). 

Some references, such as [[81, [0, use a different discrete-time model with fewer frequency 
bins, each experiencing an independent fading coefficient that repeats itself for many consec¬ 
utive symbols. We prove in Appendix |B] that the two discrete-time models are compatible. 
In Appendix jC] we show that the two models are equivalent at the continuous-time level using 
concepts of multi-carrier modulations and we provide explicit mapping of the channel coefficients 
between the two models. Therefore our results are independent of the model chosen. 

Wideband capacity is related to peakiness in the sense of the normalized fourth moment of 
the inputs, or kurtosis (bl, [[T3l . Given a stochastic sequence A{t), its kurtosis is defined as 

Eaw [|a(t)r] 


K{A{t)) ^ 


E-.,. \\n( + \\2A' 


( 6 ) 


-A(t) [|a(t)| 

where the time index (f) may be dropped if the process is stationary. By enforcing a duty cycle 
5e(0,1] on the input signal x, the system is converted into the time-alternation of an active stage 
for a fraction 6 of the time with boosted power P'=j, and an idle stage for a fraction (1—5) 
of the time. Let x be a non-peaky signal with power P and finite kurtosis k(x). We introduce 
a binary random variable cG{0, 1} to represent the use of each fading block of size TcXBc, 
where c=l means the channel block is active for signal transmission and c=0 means idle, with 
probability Pr{c=l) = 5. We reveal c to the receiver in advance, which will reduce the rate as 

C{B, 5) = I (x; y|c) = I (x, c; y) - I (c; y) = I (x; y) - I (c; y) < I (x; y), 

where 0<I (c; y)<H (c) with all equalities hold for 5=1. The duty cycle induces a new signal 

x/\/5, w.p. 5, 

0 , 


x= X 


E[cl 


with k(x) 


.EH 

Ellx|! 


E[W 

5E [Ixl 


K X 


(7) 


w.p. 1—5, 

Therefore we can effectively adjust the peakiness (in the sense of kurtosis) of signaling without 
imposing any extra constraint on the distribution of the active signal x. 


III. Bandwidth Occupancy Limit 

Our analysis is a generalization of the the SISO analysis with non-peaky signaling in [|7||. We 
extend the process to MIMO systems and to an arbitrary level of signaling peakiness through the 
tunable duty cycle parameter 5 e (0,1]. Both analyses follow four steps, represented in Fig. [2l 
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1) Find a bell-shaped lower bound R^^{SB) < C{B,S); 

2) Determine the unique maximum of R^^{6B), R^^{{6B)*); 

3) Find a bell-shaped upper bound R^^{5B) > C'(i?,(5); 

4) Determine {5B)^ and {5B)~ such that R^^{{5B)^)=R^^{{5B)~)=R^^{{5By). 

The result of [|71| shows that the capacity of a non-coherent fading channel with non-peaky 
signaling ( 5 = 1 , finite k) grows with bandwidth B only when it is below a critical bandwidth Bent, 
which falls into the range [B~, 5+]. A system operating with insufficient bandwidth B<Bcnt is 
less efficient in converting available signal energy into rate due to the sub-linear law between 
rate and SNR, and the corresponding achievable rate grows with increasing bandwidth. When 
signal power spreads over too much bandwidth B>Bcnt, the channel-uncertainty induced penalty 
grows with increasing bandwidth and the achievable rate decreases to zero as S—)-oo. Therefore, 
contrary to the wideband AWGN channel where “the deeper into the low-SNR regime, the 
better”, in the non-coherent fading channel the guideline is “enter, but not in excess, in the low- 
SNR regime”, with the optimal operation point at Scrit- Our result shows that for any B>Bcnt 
it is possible to bring the capacity back to the same optimal value, up to a small approximation 
error of order o{l/Bcnt), by imposing a duty-cycle parameter 6={6B)cnt/B and a power-boost 
P' = P/6 on the original non-peaky signaling. Moreover, in Sec. |lV]we show that this strategy 
achieves the same gap from C°° as in the peaky-signaling analysis Il8l. Il9l. 


A. Capacity Lower Bound for C{B, 5) 

As in [[7]|, our lower bound is obtained by first calculating the maximum achievable rate 
of a coherent non-peaky signaling under the average power constraint and then deducting the 
maximum rate penalty from lack of CSI. Therefore it is valid for general channel fading 
distributions and for any value of B and BcTc>Nt. The potential spatial correlation among 
different antennas is not considered here. 


Lemma 1. The achievable rate in a wideband non-coherent channel with i.i.d. fading and a 
duty cycle 5 G (0,1] is lower bounded by 


R^^{6B) = 


PN, 


1 - 


P{K-2 + Nt + Ny' 
25BNtNo 


SBNtN, 

BfTr 


log f 


p 


^ 6BNtNo 


BcTc I , 


( 8 ) 


where k = K{h) is the kurtosis of the channel fading coefficients. 








II 





Figure 2. Our four-step analysis of critical bandwidth occupancy. Substituting {SB) by B gives the original analysis by Q. 


Proof: See Appendix |Dl ■ 

The kurtosis k for many fading distributions are in the range of [1,2]. For example, as given 
in [|71, K = 2 for Rayleigh fading, k = 2—4fc^/(l+2A;)^ for Rice fading with factor A; > 0, and 
K = 1+1/m for Nakagami-m fading channels. 

Remark 1. Even though the duty-cycle constrained capacity C{B, 5) might be two-dimensional 
function of 6 and B , the lower bound is only a function of the product 6B. 


B. Maximum of 

We use the assumption that ilcT'cS>l to determine the maximum of the capacity lower bound 
. For any finite BcTc we can approximate the optimal {5B)* and the associated maximum 
rate up to an error term o ( J 


that decreases to zero as BfRc —)■ cxd. 


Lemma 2. R^^{6B) is maximized at 6B={6B)* with 


{6By = 


R^^fdB)*) > 


BrTr 


NoNt V log(5eT,) 
PN, 


(ft — 2 + Nt + Ny + o 


BcT, 


log(SeTe) / ’ 


N. 


1 - 


1 + \og{BcTc 
BrTr 


(ft-2+At+^r) logvr 


— o 


'l og(^cTe ) 

BrTr 


(9) 

( 10 ) 
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Proof: See Appendix |El ■ 

Remark 2. We would like to emphasize that the rate-maximizing bandwidth occupancy {6B)* 
is very large given the fact that the channel coherence BfTc usually ranges from a few hundreds 
to hundreds of thousands. For example, assuming 2x2 MIMO over Rayleigh fading (k= 2) with 
P/No=70 dB, we have ((5i?)*~120 MHz with capacity gap A/C°°<0.18/or i?cTc=10/ and 
(55)*~930 MHz with A/C'°°<0.03/or Bfr^=W)^. 


C. Capacity Upper Bound for C{B, 5) 

We obtain a eapaeity upper bound for the ease when ehannel is Rayleigh distributed. The 
bound, up to an error term of o{l/5B) that vanishes as 6B —)■ oo, applies to any value of B and 
BcTc>Nt and all inputs subjeet to eonstraints of average power P and signaling duty eyele 5. 


Lemma 3. The achievable rate of signaling schemes with duty cycle 5g(0,1] in a wideband 
non-coherent Rayleigh fading channel is upper bounded by 




1 - 


SBNtNo^ 


2SBNn 


PBfP, 


log(l+ 


P 


6N,BN( 


'BfPcQmm'f') 


( 11 ) 


where 0 ^ 1 , 1 = mium „ „ E is the minimum non-zero square channel gain among all 

delays and antenna pairs, and the random variable f is defined as 


/= arg minE 

i’K.n 




P 9minfK,n 

6NtNoB 



where 


2-^k=0 




( 12 ) 


Proof: See Appendix IB 


Remark 3. The auxiliary variable ip is bounded (/>0, F[f}]<l) and serves here as a placeholder 
for the minimization of the last term of the bound, which is implicitly determined by (fTB) . 


D. Critical Bandwidth Occupancy 

We obtain the range of values of 5B where the upper bound is larger than R^^{{6B)*). 
C{B,5) ean approaeh C°° within the small gap in (fTB) only if the bandwidth oeeupaney is 
eontained in an interval that grows linearly with as suggested by dB), and the error 

term o(^) in Lemma [3] can be substituted with an equivalent term 
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Lemma 4. In a wideband non-coherent Rayleigh fading channel, the maximum rate in (fTOl) is 
achievable at a critical bandwidth occupancy {6B)crit that resides in the range 


where 


{5B)- < {5BU, < ( 55 )+, 


m- 


BcTc 


No2^{Nt + N,) log Try log(S,Te) 


BcTc 




+ o 


+ o 


Proof: See Appendix O 


BfTc 


\og{BcTc 

BcT, 

logiBfTc 


( 13 ) 


( 14 ) 


E. Interpretation of the Result 

Our upper and lower bounds on C{B, <5) are all derived from the ehain rule 

I (x; y|c) = 51 (x, H; y|c=l) - 51 (H; y|x, c=l), 

where the first term eorresponds to the data transmission setup that quantifies the information 
about Hx eontained in y and the seeond term ean be interpreted as a “ehannel estimation” 
setup that quantifies the rate penalty for not knowing H. Both terms grow as 5B inereases but 
the first term grows faster when 5B is small, thus inereases I(x;y|c), until the seeond term 
“aeeelerates”. Beyond the eritieal point (5i?)crit, the seeond term grows faster than the first, thus 
erodes I(x;y|c), until the eapaeity drops to zero when 5B —)■ oo. This behavior is illustrated 
in Fig. [2] for both the eapaeity upper and lower bounds. The eoherenee time Tc and eoherenee 
bandwidth B^ of the ehannel jointly determine the relative speeds of this “raee” through their 
produet. The faetor TcBc appears in rate penalty both as the denominator outside the logarithm 
(there are Tf.Bc times fewer i.i.d. ehannel realizations than signal realizations) and as a multiplier 
of the SNR inside the logarithm (the power of TcBc signal realizations ean be eombined to 


estimate eaeh ehannel realization), leading to a eapaeity gap depending on 


Og BcTc 

BcTc 


In Fig. |3(a)| we represent the upper bound to eapaeity as a field over the 2D plane (5, B) 
with BcTc=10^ and P/Nq=20 dB. In the vertieal eut for 5=1 we also plot the lower bound 
using triangular bullets. Note that we intentionally ehoose a smaller value of P/Nq to illustrate 
the details of the transition phase in eapaeity, whieh would otherwise be diffieult to observe 
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(a) Upper bound for (S, B) and the lower bound for 5=1. 



(b) Contour plot and levels of bandwidth occupancy. 

Figure 3. Capacity upper bound over the plane {S, B) and the low bound for 5=1, with BcTc = 10^, A^t=A^r=l, and an 
intentionally chosen small value P/A^o=20dB. Range of critical bandwidth occupancy is also shown. 


with typical values of P/A^"o~70 dB [|7||. On the P-axis, we can see that for fixed values of 
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S the capacity as a function of bandwidth is bell-shaped, grows at small bandwidth, reaches a 
maximum and then decreases to zero. Fig. |3(b)| provides a better perspective on the value of 
capacity upper bounds as a function of the bandwidth occupancy SB, where the optimal (SB)* 
that maximizes the capacity lower bound and the range [(<51?) (SB)^] for the critical 
bandwidth occupancy (5i?)crit are also plotted. For different level of peakiness S, the peak 
values of capacity are the same but appear at different bandwidth B, and in fact all points with 
identical value SB have the same lower/upper bounds. Our analysis recovers previous results 
for non-peaky signals by setting <5=1, producing a finite critical bandwidth. It also recovers the 

capacity with infinite-fourth-moment signals by taking i-O, which drives the critical bandwidth 

(SB)crit 

occupancy point further into higher bandwidths satisfying lim--— = oo. 

5^0 S 

Our analysis also unveils the impact of the dimensions of the MIMO array on the non-coherent 
wideband fading channels. The maximum rate in Lemma [2l derived under the condition that 
Nt<BcTc, depends critically on channel coherence BfTc.. For example, for B^T^ = 10^ and 
B^Tc = 10®, the maximum rate (fTOl) can be approximated as, respectively. 


R(B,T,= IQ^) ~ ^(1 - Rl^/Nt+Nr), i ?( 5 , T ,= 10 ®) ~ ^(1 - O.OO50Vt+i^). 
iVo iVo 

When Nr>l is fixed, increasing the number of transmit antennas will degrade the rate, with 
the gap growing linearly with ^JNt+Nr. When Nt>l is fixed and channel is relatively flat 
(hence BcT^ is large), the rate gap is negligible for typical MIMO setups and therefore the rate 
grows almost linearly with Nr. When the channel is rather dispersive (hence BcTc is small), 
however, increasing Nr will produce a power gain that increases the rate at speed Nr but at the 
same time will bring in more channel uncertainty that increase the penalty at rate proportional to 
■\/Nt+Nr. Therefore using too many receive antennas will hurt the achievable rate. For example, 
the maximum rate peaks around iVr=40 for BcTc = 10^. It must be noted that our analysis is 
accurate for conventional MIMO systems with Nt < BcTc, and extension to high-dimensional 
MIMO is out of the scope of the current paper. 


IV. Unified Capacity for Peaky and non-Peaky Signals 

In this section we will show that the peak rate R^^((SB)*) in (fTOl) . which is derived by 
combining non-peaky signaling analysis |I3 and tunable peakiness through duty cycle (5 g( 0, 1], 
approaches within the same gap as in the unconstrained capacity C(B) analysis using a 
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generalized polynomial rate approximation Q obtained via peaky signaling analysis BUl. 
We first replaee 5, a free parameter in our model representing the duty eyele, with an assigned 
value 5=SNR^“" as in where SNR = P/{BNq) and a is the exponent that determines the 
wideband slope. This substitution will show that, when the channel coherence length is large, i.e., 
BcTc 1, the gaps to C°° in (fTOl) and in Q have the same value at points {B, 5 = SNR^“"). 
Furthermore, we show that the sufficient and necessary conditions on the coherence length BcTc 
to approach C°°, proved in lIH Th. 1-Th. 3], can also be established using our results. Once we 
have established that the results are equivalent, the opposite path can be taken and use the values 
of a obtained in @ to calculate a new range of the critical bandwidth occupancy in closed-form 
expressions. We discuss the relationship of the two expressions, which have minor differences 
in the error terms of the calculation of a, reveal a trade-off between accuracy and resolution in 
Il9l , and demonstrate that the two methods represent the same optimal rate. 


A. Different Analyses Show the Same Results 

The analysis in IQ obtains a necessary and sufficient condition on the coherence length of the 
channel, BcTc, to guarantee that capacity is above a polynomial of S'NR=j^ as B^oo with 
specified peakyness 5 = SNR"“^. This result is given in |l9l Th. 3], which is rewritten in the 
next lemma for easy reference. The result is valid for arbitrary BcT^, but the necessary condition 
to approximate is akin to requiring that BfTc be large. 


Lemma 5 (Th. 3 JH). For any aG(0,1] and eG(0,a) the capacity of a Rayleigh block-fading 
MIMO channel with coherence time Tc, coherence bandwidth Be, and average signal to noise 
ratio SNR=-g^ is 

> TVrSNR - + 0(SNR^+“+"), (15) 


if and only if there exists a cTG(0,e) such that 


\Ne + N,y 


N‘^ 

* SNR-2('^+“). 


(16) 


Recall that in Sec. Unl our rate lower bound in ([8]) contains three terms, the wideband capacity 
C°°, a non-linear rate penalty due to log(l-l-SNR), and a rate penalty due to lack of CSIR. Below 
the optimal bandwidth occupancy (SB)*, the third term of ([8]) is smaller in absolute value than 
the second. Replacing the third term by the second term and substituting (5=SNR^“", q;G(0, 1) 
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into dl]) produces the following suffieient eondition in terms of the bandwidth oecupaney 5B, as 
stated in Corollary [TJ 


Corollary 1. If 5B < {6B)*, the achievable rate is lower bounded by 


C{B,5) > 


PN, 


1 - 


P 


(k— 2+iV-t+A^j.) 


(17) 


On the other hand, above {6B)cnt, the third term of ([8]) is greater than the seeond. This means 
that C{B, 5) is smaller than (fTVl) . whieh leads to the necessary eondition in Corollary [2l 


Corollary 2. In Rayleigh fading (k= 2), if 6=SNR^ " and if 


C{B,6) > 


PN, 


1 - 


p 


(iVt + iVr 

N, 


(18) 


then the bandwidth occupancy satisfies 6B < {SB)~^. 


Now we ean use the neeessary eondition in Corollary |2l and the suffieient eondition in 
Corrollary [T] on bandwidth oeeupaney 5B to prove the suffieient and neeessary eondition (fT^ . 

Proposition 1. Corollary [7] implies the sufficient condition (fT^ for Lemma |5l 

Proof: Substituting 5=SNR^““ and k=2 into dH), we ean rewrite 6B<{6B)* in Corollary [U 


as 


BcTr, > 




SNR-^“(iV, + iVt) log(ReTe). 


(iVr + iVt)2 

iSnee {N^ + Nf \og{BcTc) is a eonstant and R —)■ oo, we have {N^ + Nt) \og{BcTc) < SKR”^*^ 
for any e > 0. Therefore, it is also suffieient to have 

N? 


BrTr > 


-SNR-2(“+^), 


{K + 

whieh is a suffieient eondition that Lemma [5] transforms in the upper limit of a < e. 
Proposition 2. Corollary \2\ implies the necessary condition dlH) for Lemma |5] 

Proof: The neeessary eondition 5B < in Corollary [2] can be rewritten as 


> (iVr + iVt)2 




dlogTT 
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Therefore we ean express the neeessary eondition that Lemma [5] sets as lower limit of cr > 0, 

-SNR-2", 


BT , 


as long as log(RcT'c) > i-e., whieh is always satisfied in wideband 

fading ehannels where BcTc is very large, and thus BcTc > vr^. ■ 


Remark 4. From Proposition [7] and Proposition |2] it is not surprising that the power gain term 
\og{BcTc) was lost in because this sub-polynomial variation of the result has been “buried” 
in the range of valid exponents e of the error term 


B. New Bounds on ((5R)crit using the Subquadratic Polynomial Rate Approximation 

In our analysis, Lemma [2] preseribes a near-linear-in-power eapaeity lower bound whieh ean 
be achieved by all signaling schemes with (5, B) as long as the bandwidth occupancy 5B equals 
some constant (5i?)crit- Our analysis does not provide the exact value of {5B)cr\t, but rather 
bounds it within [{6B)~, {6B)~^] in Lemma HI On the other hand, the result in |I3 Th. 3], 
reproduced here as Lemma [51 prescribes an entire family of parametrized bounds where the 
parameter e controls both the error term of the generalized Taylor expansion and the resolution 
of bounding brackets around {6B)cnt- Corollary [3] makes this explicit. 


Corollary 3. The necessary and sufficient condition (ITbl) of Lemma |5| shows that coherent 
capacity C°° is approached by transmitting signals with bandwidth occupancy SB within the 
limits 




SB > 


No N, 

P fN, + Nt 


No 


Nt 

Proof: Substituting BfT^, 5=SNR^-“ and SNR= 
by the fact that (j>0 and a<e, respectively. 


^ = (SB) 


mm 

e 


(19) 

( 20 ) 


BNo 


into ([T6l) . we can obtain ([T9l) and (l20l) 


Therefore for a given a, which controls the level of peakiness S and determines the wideband 
slope, we can observe a clear tradeoff, parametrized by ee(0,Q!), between the accuracy of the 
Taylor polynomial and the resolution of the bandwidth brackets: 
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1) The accuracy of the eapaeity lower bound calculated in Lemma [5] is determined by the 

ratio between and the error term The larger e, the better the 

approximation, since the error term will vanish faster as R —)■ cxd. 

2) The resolution of the interval where (^R)crit is contained, [(55)“™, is determined 

by the width of the interval. The smaller e, the better the resolution, as the lower boundary 
((5R)™™ will increase and become tighter when e becomes smaller. 


C. Comparison of Critical Bandwidth Occupancy Estimators 

So far we have characterized the critical bandwidth occupancy ((5R)crit in two different 
ranges: by a pair of bracket (5R)+], in our analysis in Lemma |4l and by a parametric 

interval [((5R)“™, (55)™'^^], derived from |I3 Th. 3] in Corollary [3l To explore the relationship 
between the two estimators, we compare the difference in the estimated value of a that each 
analysis produces. We do this because the exponent a provides a unique relation between B and 
5=SNR^“"=(j^)^“", allowing for scalar comparison of the methods. 

We begin by representing a according to dH Th. 3]. From (fT^ in Lemma [51 for given values 
of the coherence block length BfTc and bandwidth B G [(5R)™™, can be written as 


log( 


a 


a = 


TV? 


BfTf) 


21og(SNR-^) 


( 21 ) 


From the fact that a > 0 we get 


log( 


(TVt+Af,)^ 

TV? 


BfFf 


21og(SNR 


-It 


( 22 ) 


and from the fact that a < e < a we get 


Q^max 

a > max —-—, 


e) , where 


fe)^ 


log( 


{Nt+N,f 

TV? 


RcTe) 


e. 


(23) 


21og(SNR-i) 

Note that when e decreases, amin(e) increases such that the range of a becomes smaller but at 
the same time the error term 0(SNR^’''"’''^) vanishes more slowly: improving the resolution of 
the bandwidth occupancy range comes at the price of decreasing the accuracy of the capacity 
polynomial approximation. We can make an approximate selection of e such that polynomial 
error term is in the order of a p-percent of the term SKR^"*"", i.e., finding e(p) such that 

100 , 


SNR 


l-\-OL 


> 


P 


-SNR^+“+^(^’). 
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This generates a family of narrower estimated margins [amin(p), «max] parametrized by the pre- 
seleeted error pereentage p% by raising the lower braeket. 

On the other hand, we ean bound a using the eritieal bandwidth oeeupaney interval in Lemma |4] 
in eombination with 5=SNR^“". With 6B={SB)~^ we get 


log( 


4 (Aft + Afr ) log TT BcTc 


a = 




\og{BcTc )' 




log( 


(Nt+N,)\og{B^n) s 
4 log TT ^ 


21og(SNR-^) 
and with 6B = {SB)~ we get 


21og(SNR 




21og(SNR-^) 


(24) 


log( 


BnTr 


a = 


4(Art+Af,)log7rlog(Bcrc)' 


log(^^^^±^R,T,) log(4 log TT log(ReTc)) 


2 log(SNR- 


2 log(SNR 




21og(SNR 


-1^ 


(25) 


Reeall that for any e > 0 we have (iVt + N^) log(i?c^c) < SNR This means that we 
ean show that ctmax >«'*'> amin(e), and the interval between the three vanishes as e—)-0. 


Remark 5. All the results coincide in that a<x.\og{B(.Tc), making capacity of channels with low 
BfTc approach their wideband limit very slowly with SNR^O and channels with high BfTc 
approach the wideband limit faster. This is the main intuition of the results in non-coherent 
channels approach the coherent channel capacity when coherence length is large enough. 


D. Illustration 

We plot the eapaeity lower bound on the plane (5, B) in Fig. |4] for BcTc=lif (first graph) 
and for BcTc=l^^ (seeond graph). The peak eapaeity is aehievable in a region with eonstant 
produet 5B, starting at relatively large bandwidths, and both estimations of the optimal region 
are narrow. The ehoiee of e determines the polynomial lower bound and therefore the range 
[amin(e), amax]- We ean generate a set of estimations Q!min(e) by fine-tuning e within the range 
(0,a), as shown by the eurves eorresponding to amin(e) with e=a/2,a/4, respeetively. Note 
that the eonservative ehoiee e=a leads to the widest possible range for [Q;max/2, ttmax]- On the 
other hand, the resolution of the estimators from our own analysis [a“, a^] depends only on the 
value of BcTc, and its range beeomes smaller as BfTc inereases. 

Sinee the resolution of the estimation by [a~^a^] relies on BcTc and the relative margin of 
[ttmin, ttmax] depends on e, we show in Fig. [5] the evolution of the two boundary methods with e 
and BcT- The method IIH produees the highest upper bound Omax that does not ehange, and a 
family of lower bounds Q!min(p) depleted in the figure for errors of 1% and 10% and its lowest 
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B T =10® 



Figure 4. Capacity lower bound on the plane {S, B) with P/A^o=20dB, _Bc 1^=10® (first graph) and BcTc=10* (second 
graph). Curves with Q:min(e) are generated with e=a, a/2, a/4, respectively. 


bound Note that at low coherence length, the limit ctmin > amax/2 makes it 

impossible to select values of e corresponding with a polynomial accuracy of 1%, and then 10%. 
This shows that the polynomial rate with peaky signaling in dUl also displays a gap from C°° 
decreasing with B^Tc. On the other hand, the critical bandwidth occupancy method produces 
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Sublinear exponent a bounds versus B T SNR=-20.0dB 

c c 



Figure 5. Evolution of y ) versus B^Tc with SNR = —20dB. 


boundaries that are loose at low coherenee length but improve signifieantly when this parameter 
grows and that do not pay for tightness a priee in aeeuracy of the polynomial approximation. 

V. Conclusions 

In this paper we have unified the study of the rate approximations to for peaky and non- 
peaky signaling in non-coherent wideband fading channels where energy rather than spectrum 
is the limiting resource. We have generalized the critical bandwidth analysis 0 to families 
of signaling schemes with varying bandwidth B and transmission duty-cycle 5 g( 0, 1] to allow 
arbitrary levels of signal peakiness. We introduce the metric of bandwidth occupancy to measure 
the average bandwidth usage over time and define it as 5B, the product between the bandwidth 
and the fraction of time it is in use. Our main result shows the existence of a fundamental 
limit on the bandwidth occupancy in non-coherent channels for any level of signal peakiness. 
For all signaling schemes with the same bandwidth occupancy, as the bandwidth occupancy 
approaches its critical value ((5i?)crit7 rates converge with the same asymptotic behavior to the 
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same almost-linear in power value (measured in nats/s) 


C{B,6 = 


(SB) 


crit' 


B 


> 


PN, 


1 + \og{BcTc 
BrTn 


■(k - 2 + iVt + N^) logTT 


where Tc is the eoherenee time and Be is the eoherenee bandwidth. The rates deerease to zero 
as the bandwidth oeeupaney goes to infinity. Moreover, we provide upper and lower bounds to 
this eritieal value. The bounds have the same growth with B^Tc and and they only differ on 
a eonstant term. 

To eharaeterize the relation between the eapaeity with a tunable peakiness eonstraint C{B,5) 
and the uneonstrained non-eoherent wideband eapaeity C{B), we rewrite the above eapaeity 
expression as a polynomial equivalent to the analysis in [|9]|. We have reeovered the results 
in Th. 1-Th. 3] and obtained the almost-linear polynomial expressions for eapaeity in the 
limit 5B —)■ ((5i?)crit with a dominant sub-linear term SNR“. As the bandwidth oeeupaney 
approaehes the limit, eapaeity approaehes the power-limited wideband limit with a speed of 
eonvergenee determined by whieh approaehes that of eoherent ehannels as BcTc —)■ oo. 

The fundamental nature of the bandwidth oeeupaney measure refleets the faet that eapaeity of 
any signaling seheme is eontained within the same bounds as long as the produet SB is eonstant. 

Within this framework, limited bandwidth transmission with non-peaky signaling and unlimited 
bandwidth transmission with peaky signaling, whieh have been treated as very different sehemes, 
are shown to be merely two extreme points in a eontinuous range of transmission strategies 
within the same bounds as long as they have the same amount of bandwidth oeeupaney. This 
suggest that for the praetieal goal of operating at a rate very close to (7°°, all pairs (R, S) 
with the optimal oeeupaney do not exhibit signifieant differenees. Aehieving eapaeity, i.e. the 
supremum rate, may on the other hand only be possible in some speeifie distributions. The 
seleeted peakiness S = SNR^““ in [[91 beeomes invalid if SNR>1 (as ^ < 1 by design), whereas 
our model determines peakiness through 5R<(5R)crit, a quantity that is well defined for all 
values of SNR. This gives the intuition that below the eritieal point it would be questionable 
to elaim that the frequeney-seleetive ehannel is in the wideband regime, and therefore regular 
non-peaky transmissions with full bandwidth oeeupaney must be employed. Beyond the eritieal 
point, both signaling sehemes provide the same eapaeity limit. 

We have shown that most of the advantage of peaky signaling stems from harnessing power for 
long periods of time to transmit some infrequent flashes with boosted power, without eneoding 
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information in the position of the aetive symbols as in ON/OFF modulations. Moreover, this 
power boost does not in faet outperform non-peaky transmission with the optimal bandwidth, 
whieh means that in practieal systems the amount of peakiness and the bandwidth may be ehosen 
at will as long as the maximum oeeupancy level is respeeted. 

Our analysis has some limitations. Firstly, the potential spatial eorrelation among MIMO 
antennas is not accounted for. Secondly, although our capacity lower bounds are valid for 
general fading channels, our upper bound and critical bandwidth occupancy expressions assume 
Rayleigh fading. Besides, the performance of a signaling system with practical channel estimation 
techniques ll22l . peak constrained signals ifTTll . [fTSl . lfT9ll . Il2^ . finite modulation options, and 
non-ideal decoders may be degraded as compared to the theoretical bounds provided in this 
paper. 


Appendix A 

Justification of Our Fading Modef Choice 


As a general case, a wireless channel is modeled as a set C of paths, where each path £ G £ 
is defined by a group delay ti, a phase of arrival 9(,, and an impulse response hi{t). For a pair 
of antennas (m,p) with received signal and transmitted signal we have 




(t) = (t) * ^ {t - (t) = (t) * (t) + (f), (26) 


e&c 

where is the AWGN noise, and the channel delay spread D and coherence time are 

determined by the aggregate channel impulse response Traditionally, are scalar 

gains or narrow pulses that can be approximated by the Dirac delta function, in which case the 
set C would be a sort of “ray tracing” of perfect reflections of the signal with a scalar gain. 
However, recent mm Wave meassurements have found much higher delay-spread values ll2^ than 
those predicted in ray-tracing calculations ll25l . This may be due to rich scattering from small 
objects in mm Wave fading channels, which are not so sparse in practice. This is due to the fact 
that, although there are few arrival direction “clusters”, in each cluster energy arrivals spread 
along many angular directions [l26l . Therefore each arrival direction sees the additive effect of 
a large number of scattered reflections, not a single path, and each has a delay spread, 

instead of a scalar channel gain. The construction of discrete-time system models falls into the 
following three regimes depending on the sampling rate: 
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• Sampling at low rate, all the energy in the delay spread D would be captured by 

a single sampling interval, so the resulting discrete channel would be a scalar coefficient, 
which is approximately Gaussian distributed due to the law of large numbers. This is called 

the narrowband, or frequency-flat channel. 

• Sampling at higher rate would make the energy in D be captured in multiple sampling 
intervals, each with an independent scalar coefficient. This is called the wideband channel, 
or frequency selective with rich scattering environment. 

• The third regime occurs when the number of sampling bins is much larger than the number 
of paths in C. The sampled channel coefficients are sparse and not Gaussian distributed. 
This is called the ultra-wideband. 

We consider the wideband fading model to be relevant in mm Wave communications where 
rich scattering and longer delay spread was observed lf24l . Our discrete equivalent channel is 
derived from the propagation described above by employing the classic framework of a Nyquist 
sampling at frequency B, the consideration of frequency-domain signaling with a iT-point DFT, 
satisfying K = BTc, and a cyclic prefix of negligible duration M = B/B^ = K/BcT^ -C K. 
The same channel model is employed in 0. 

Appendix B 

Compatibility With Another Common Model 

In dUl the signals are divided into a set of M = narrowband channels (a.k.a. frequency 
bins) with encoding symbols defined with a symbol period of 1 / B^. Each narrowband channel 
can be perfectly sampled at a rate of just 1 sample per symbol period, and there are M parallel 
frequency bands that produce M samples per symbol period. In this scheme multiple symbols 
see the same channel realization and the channel coherence length is a block of Lc = TcBc 
consecutive symbols. By indexing with m the independent frequency bins and with I the 
consecutive periods on the same channel block realization, we get the model 

y[m,f] = H[m, f]x[m,f] -l-z[m,f], (27) 

where H[m, £] remains unchanged for £ = Lc. To exploit channel coherence, the encoding 
process must design the transmitted signal for the Lc consecutive symbols jointly, and the 
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encoding model is represented with matrices as 

Y[m] = H[m]X[m] + Z[m], (28) 

where the dimensions are Lc = x Ni){Ni x Lc). 

In this model, for every encoding interval of length Lc and across all M frequency bins there 
are a total of MLc=K complex valued coefficients. Therefore, this channel model provides 
exactly the same number of signaling dimensions for transmission as the model we have derived. 
But the representations of the channel variation are different. In this model there are fewer channel 
coefficients, each of them is i.i.d. and identically repeated for every Lc consecutive symbols. 
Whereas our derived model dU) supports any type of channel correlation, not only repetition, 
as long as there are K correlated coefficients generated by a fraction 1/BcTc of independent 
random variables. It is possible to represent the system model (1281) with repeated identical channel 
coefficients in our derived model format by replacing the matrix notation H[m]X[m] with our 
vectorized notation Hx where H is a block-diagonal matrix with the values of H[m] in its main 
diagonal and zeros in the upper and lower triangles as in dH). 

Appendix C 

Equivalence in Signaling Representation 

Our channel model uses Nyquist sampling at the full B and therefore it is able to represent 
any signal with this bandwidth without loss. For the sake of completeness we will propose the 
exact formulation to implement a valid signal in the model of dH (hereafter, filter-bank model) 
with our model (hereafter, OFDM model) using only preprocessing linear matrices. With this 
we show that any signal possible in the filter-bank model can be transmitted through the OFDM 
model, and therefore capacity results in our model are fully compatible. 

Without loss of generality, let us assume a SISO channel and unit power to simplify notation. 
Assume also that the integers K=\TcB^, M=\B/Bc^ and Lc=\TcBc^ are satisfied exactly so 
we may use simply K=MLc. In continuous time, the filter-bank model is represented in 

oo /Lc-1 /M-1 \ \ 

>■(*)= X 5{t-iTc) * I ^ 5{t-i/Bc) * I ^ hi[m, i]xi[m, i]smc{tBc)e \ \ 

i=—oo V £=0 \m=0 / / 

(29) 

where multiplication by hi[m, ^] and Xi[m, ^] assigns the scalar value received in each frequency 
bin m e {0,..., M—1} and in each transmit symbol period £ G {0,..., Lc—1}. 
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We separate the eneoding for eaeh ehannel bloek realization indexed by i, drop the index, and 
use the faet that the ehannel eoeffieient in eaeh frequeney bin remains the same for all symbols 
to take away the index I from h[m\. This gives 


Lr— 1 


^M-l 


y{t) = ^2 * I ^ /7[m]a:[m, f']sinc(f-Bc)e-^^™^‘'* I +z{t). 


(30) 


£=0 


^m=0 




With this eontinuous-time signal, we apply Nyquist sampling at rate B to generate K=BTc sam¬ 
ples per sequenee. Notiee that for integer M=B/Bc, the diserete sine funetion is sinc[n/M] = 
sinc(^^) and the delta delay on index n is iM. We ean represent (l30l) by 


Lc — 1 


^M-1 


y[n\= ^ 6[n—iM] * ^ h[m]x[7n,i]smc\ — \e-'"“i< 


,j2TT^mLc 


e=o 


vm=0 




+z[n], n = 0,..., 1. (31) 


We eompute the if-point DFT, 


= -i2^T7 


i-<^h[rn\x[rn^l]iec[k/Lf. — m\+z[k]. 


(32) 


i^m 

The reetangular window equals one only when [k/Lc\=m. By representing k<K as k = 

A I 7 / r I 1 A 


u* Lc + V with u = \_k/Lc\ and v — k mod Lc, we obtain 


Lc-l 


h[u] ^ x[u, i] + z[k], with 
e=o 


u=[k/Lc\, 
v=k mod Lc 


Now we ean see that the sum is aetually the nth element in the Lc-IDFT of the sequenee x[u, £]. 
Sinee the IDFT of a sequenee a=(ai... a^^)^ ean be written as a matrix produet IDFT(a) = Fa, 
We ean represent the system model the same way as our matrix ehannel notation as 


y = H<Fx -f z, 


(33) 


where H for SISO is a if x if diagonal matrix with its fc-th diagonal element /i[n], x is if x 1 
with = x[u,v]. The Lc-IDFT is eomputed by the bloek-diagonal square matrix 

\ 



^ F 


0 

$ = 





1 " 


F 


(34) 




This shows that any ehannel of the filter-bank model ean be represented by the OFDM model 
using a ehannel matrix H = H<F. The reeiproeal eompatibility ean be proven by taking a 
preeoding DFT matrix at the transmitter, x = $lx, whieh leads to 

y = Hk = = H$$+x = Hx. 
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The multiplication by <1)1 is unitary, so if the OFDM model uses x ~ CM(fx, S) and $ is a 
full-7F-rank square orthonormal matrix, then <l)lx ~ CAA(<l)l/x, <1)1S<1)). Gaussian distribution is 
maintained when the channel model is changed, and the mutual information results for both 
channel models supported by our bounds based on Gaussian inputs are completely equivalent. 


Appendix D 
Proof of Lemma 1 

Since the receiver knows which phase the duty cycle is in (e.g., scheduled according to a 
pseudo-random sequence), the rate can be determined via the chain rule 

(x; y|c) = (x; y|c=l) + (1-5) ■ 0 = (x, H; y|c=l) - (H; y|x, c=l), (35) 

J- r J- r J- r J- r 


where the first step comes from the fact that Hx=0 in the idle block (c=0) and Pr{c=l) = 5. 

p ' 

SBNo ■ 


During the active block the input follows a Gaussian distribution CM{0, MM) and the first term 


in (135]) can be lower bounded by 

(x, H; y |c=l) > (x; y |H, c=l) = 5 x Eh 


Tc 6BNtNo 


HHl) 


(36) 


where the first step is from the non-negativity of mutual information, and the second is due to 
independence of channel coefficient H in each subcarrier and transmit antenna. Furthermore, 
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>SBEu 


(pPM 
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{^PN, 
~ No 

PNr 


Ptr(HHt) / P ytr(^(HHl)= 
6BNtNo ~\6BNtNo' 


P/m 

2MN/No 




\ht,M\hu,r 

t,r t^u,r 




t^r^v 


ty^u,r^v 
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0 


P {MMk + NtM{M-l) + MM{M-1)) 
26BMNiNo 

(k — 2 P iVt P M) 


26BNtNo 


(37) 


where Aj are eigenvalues of HHl with il=[ht,r]NrxNt representing the diagonal blocks of H, 
and Mr is the (r, f)-th element in H. Equation (a) comes from the fact that H[fc] are identically 
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distributed for all A;=0,... ,K—1, (b) is due to log(l+a;)>a;—a;^/2 for x—)-0 and the fact that 
Ai=tr(HHl) and A^=tr((HHl)^). Equation (c) is by careful rearrangement. Equation (d) 
comes from Eh [|/ip]=l, Eh [|/i|^]=k, Eh [h]=0, and independence of matrix entries. 

To upper bound the second term we choose H to be Rayleigh fading (with the maximum 
entropy) and interpret x as a pilot signal that gives side information between H and y. 

I (H, y|x, C 1) ^ I (HGaussian; y|^Pilots Signal) C l)^I(H;y|x,c=l). (38) 

An example for channel estimation would be a system where the pilot signal transmitted on 
antenna m is a uM times delayed version of the signal on antenna 1. After transmitting K pilot 
symbols, at each receive antenna a A'-equation MAt-unknowns linear estimation problem is 
established and can be solved using the MMSE estimator. 

Eet be the MNt X MNt diagonal matrix containing in its uM+m diagonal element 
QuM+m = E (the gain of the m-th channel tap in the {u,v) transmit and receive 

antenna pair), and let S be a A x MNt circulant matrix (MNt<K) containing mod x 

in its (i, j)-th coefficient, where is unit-power pilot signal. Notice that the mention 

of pilot signals here is to upper bound a mutual information term, rather than implementing a 
practical channel estimation as required in a coherent receiver. Exploiting the fact that channel 
estimation is carried out on each receive antenna concurrently based on the hypothetical pilot 
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signal E from all transmit antennas, we get that the upper bound results in 


Nr 


6 7 6 V-^ 
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V MNtNo' 
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log I 1 H- 

^ V 6BNoN, 


{BtTc 


(39) 


where (a) stems from the AM-GM inequality and that ehannel gains between all antenna pairs 
are i.i.d, (6) is due to the faet that H is a eireulant matrix, whieh has the same eoeffieients shifted 
aeross all its eolumns, so its eigenvalues are the DFT eoeffieients of the eolumns, (c) is Jensen’s 
inequality, (d) derives from the faet that x has unit power, and (e) is due to the upper bound of 
squared ehannel eoeffieients J2n=i 9n ^ and (/) uses ^ = BcTc. 


Appendix E 
Proof of Eemma 2 

Taking partial derivative of ([8]) w.r.t. the produet 6B, we obtain 


dN^{6B) PNr 

d{6B) 


P{N-2+Nt+Nr) NoNt 
2{6ByNtNo PB,n 


log ( 
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V {6B)NtNo 
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SB (IH_ pMELN\ 

on 1 I 


(40) 


Near the maximum of N^{5B) the term -jj^^BcTc is either 3>1 or ~1; beeause N^{5B) is 
already approaehing zero if BcT^ <C 1. This means we ean make the approximation 


{SB)No - 

{k — 2 + iV-t + Nr) 


log(l + 


(41) 


2 {P/{NtNo{5B)*)f B.N 

whieh solves as dH). Evaluating N^{6B*) and using the same inequality in [|71 produees (fTOl) 
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Appendix F 
Proof of Lemma 3 


We upper bound the first term in (1351) enforeing signal bandwidth B and duty eyele 6. 

^I(x,H;y|c=l) = ^h (y|c=l)-^h (y|H, x, c=l) 

= (Hx + z|c=l)-;^h (z) 

S f , P \ (5 

< -h (CAf(0, Ij + BNo)] --h (z) (42) 

where (a) is from the definition of mutual information; (b) is from the ehannel model; (e) 
eomes from the faet that z is independent of x and H, and h (Hx + z|c=l) is maximized by a 
Gaussian distribution under the power eonstraint j + BNq. Use the approximation log(l+a:) = 
x—x^/2 + o{x^), we ean rewrite (l42l) as 


5 PN 

—I(x,H;y|c=l) < 

J-C jVo 


1 - 




(43) 


25BNo_ 

For the seeond term of (1351) . with the Rayleigh fading assumption, the inequality in (IMI) is 
met with equality. From there on, upper bounds are found by taking a eouple of minimums in 
the argument of the logarithm. 


^ I(H;y|x,c=l) 
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> -E 
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(44) 


Equation (a) is due to 5 fmin= E is the minimum element in the diagonals 

of and among all p’s, and (6) stems from the relation between determinant and eigenvalues. 
Sinee H is a KxMNt eireulant matrix eontaining the power normalized veetor x/\/P in its first 
eolumn, the n-th eigenvalue of is given by 



spK-l 2 ^ 

2^k=o \ 
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P^K,ni 


n = 1,..., MNt. 
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Since E [iiK,n\ < 7 ^ I Yl,k=Q ^^1^ — ^ owing to the power eonstraint E [x] < P, we obtain (c) 
by the faet that B^Tc < K and by the definition of -0 in (12). Moreover, we have -0 > 0 beeause 
the rate penalty of non-peaky inputs in aetive eyeles is non-zero (51 (H;y|x/\/5) > 0). 


Appendix G 
Proof of Eemma 4 


We define sueh that 




(SB)±No 


BrTr 


BcTc 


Substituting (l45l) into (fTTl) we obtain that 
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We separate the logarithm in two parts 
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Sinee ^[-0] < 1, the third negative part is also o{ 
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We will make this upper bound equal the aehievable value in (fTOl) . whieh leads to 
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By making ehange of variable T = f2/At we get 
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With K = 2 for Rayleigh fading, +1) > 1. We obtain the following two roots of (l50l) 



(51) 


It is ready to see that 



(52) 


Substituting them baek in (1451) we get the points (SB) and (SB)^ as shown in (fT4l) . Therefore 
the true aehievement of the maximum ean only oeeur in the range {5B)cnt £ [(55)-, (55)+]. 
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